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ABSTRACT 

When binary black holes are embedded in a gaseous environment, a rotat- 
ing disk surrounding them, the so-called circumbinary disk, will be formed. The 
binary exerts a gravitational torque on the circumbinary disk and thereby the or- 
bital angular momentum is transferred to it, while the angular momentum of the 
circumbinary disk is transferred to the binary through the mass accretion. The bi- 
nary undergoes an orbital decay due to both the gravitational wave emission and 
the binary-disk interaction. This causes the phase evolution of the gravitational 
wave signal. The precise measurement of the gravitational wave phase thus may 
provide information regarding the circumbinary disk. In this paper, we assess the 
detectability of the signature of the binary-disk interaction using the future space- 
borne gravitational wave detectors such as DECIGO and BBO by the standard 
matched filtering analysis. We find that the effect of the circumbinary disk around 
binary black holes in the mass range QMq < M < 3 x IO^Mq is detectable at a 
statistically significant level in five year observation, provided that gas accretes 
onto the binary at a rate greater than M ~ 1.4 x 10^'' [g s"-^] ^"^(M/IOMq)^^/^^ 
with 10% mass-to-energy conversion efficiency, where j represents the efficiency 
of the angular momentum transfer from the binary to the circumbinary disk. We 
show that O(O.l) coalescence events are expected to occur in sufficiently dense 
molecular clouds in five year observation. We also point out that the circumbi- 
nary disk is detectable, even if its mass at around the inner edge is by over 10 
orders of magnitude less than the binary mass. 

Subject headings: black hole physics - accretion, accretion disks - gravitational waves 
- stars: evolution 
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Introduction 



It is widely accepted that coalescing binary compact objects, i.e., those consisting of 
black holes, neutron stars, or white dwarfs, are important sources of gravitational waves 
(GWs). Ground-based detectors with higher sensitivity are currently under construction. 



e.g., Advanced-LIGO, Advanced- VIRGO, and LCGT , whi 



space-born e interferomete rs such as LISA (jPanzmann 



e there are plans for t 



ne 



and BBO ( jPhinney 



19971), DECIGO fISeto et al. 



20011 ) 



20031 ) ■ The significant improvement of sensitivity will make it possible 
to directly detect GWs in near future. LISA has its best sensitivity at around 10~^ — 10~^ Hz 
and one of its main targets is a merger event of binary massive black holes. On the other 
hand, DECIGO/BBO has its best sensitivity at around 0.1 — 1 Hz. Its ultimate goal is 
to detect primordial GW backgrounds, but it has promising astrophysical sources such 
as binary stellar-mass black holes. It is also sensitive to merging signals from binary 
intermediate-mass black holes. 

Recently, much attention has been paid to the surrounding environment of the 
GW sources. Several types of possible electromagnetic (EM) signatures associated wit h 



GW emissi ons have been proposed, such as afte rglows flMilosavljevic fc 



precursors (IChang et a 



2007 



Bogdanovic et al. 



2009 



2008 



Bode et al. 



Phinnev 



20101) and periodic emission (iHavasaki et al. 



MacFadyen &: Milosavljevic 



2008 



2005) 



Cuadra et al. 



20091 ) in 



the context of a massive black hole coalescence. Those EM signatures will be also very 
useful in localizing the position of a GW source in the sky and in identifying its redshift. 
By combining the redshift information with a luminosity distance determined by the 
GW signal, one can obtain a distance- re dshift relation. This will be one o f the best 



observational probes of the dark energy (jSchutz 



1986 



Holz fc Hughes 



20051). However, 



there are not so many works discussing EM signals from coalescing binary black holes in 
the stellar /intermediate mass range 6 — IO^Mq. Obviously, no black holes can emit EM 
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radiation, unless they interact with their environmental gas. Thus, our interest is naturally 
led to binary stellar /interstellar-mass black holes in a dense gaseous medium. 



Narayaru (120001 ) investigated the effect of hydrodynamic interactions between an 
advection-dominated accretion flow onto a supermassive black hole and inspiraling stars on 
GW signals. This effect was shown to be negligible as long as the advection-dominated 
accretion flow is concerned. The effects of the interaction with a gaseous disk on 
the GW phase evolution have been also discussed in the context of an extreme mass 
ratio inspiral (EMRI), whe r e a compact obi ec t inspirals into a sup ermassive black hole 



(IBarausse fc Rezzolla 



2008 



Kocsis et al 



2011 



Yunesetal 



201 Ibl ). In this case the 



situation is more similar to a planet migration. It has been claimed that the existence of an 
accretion disk around a supermassive black hole is detectable for a certain parameter range 
with LISA. 

In contrast, a binary composed of black holes with nearly equal masses residing in 
a gaseous environment will form a triple disk system: two accretion disks around the 



respective 



2007 



Dlack holes and one circumbinary disk surrounding both of them ( iHayasaki et al. 
20081 ). Then, the binary black holes and the circumbinary disk should mutually 
interact, which inevitably affects the orbital evolution of the binary. The gas surrounding 
the binary acts on the orbital motion through a drag force (tidal/resonant interaction), 
leading to the angular momentum transfer from the binary to the gas. Conversely, the 
angular momentum is carried from the circumbinary disk to the binary by gas accretion. 
Since DECIGO/BBO will detect many cycles of GWs from the binary inspiral, a rather 
tiny correction in the GW phase is detectable. Therefore signature of the circumbinary 
disk can arise in the GW signal even if the circumbinary disk is extremely less massive by 
many orders of magnitude than the binary. In this paper, we give an estimate how large 
the mass accretion rates of the circumbinary disk should be, so that we can detect the 
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effect of the binary-disk interaction on the GW phase evolution at a statistically significant 
level with DECIGO/BBO. Note that we evaluate the effect of the binary-disk interaction 
at the first order in a mass accretion rate of the circumbinary disk throughout this paper. 
Although our current understanding of the binary-disk interaction is poor, we assume here 
that we knew how the GW phase evolution is modified and hence the standard matched 
filtering analysis is applicable. Based on our estimate of the mass accretion rate required 
for detection, we also give a prediction for the number of merger events under detectable 
influence of the circumbinary disk. 



2. Measuring mass accretion rates with GW detectors 

In this section, we first briefly explain the basics about how the phase of GWs from 
a binary evolves as a result of the radiation reaction due to the GW emission. Second, 
we present the picture that we have in mind about the interaction between the binary 
and the circumbinary disk, providing a quantitative estimate for the correction due to 
the binary-disk interaction in the GW phase evolution. We relate the magnitude of the 
correction to the mass accretion rate of the circumbinary disk. Then, we evaluate the 
determination accuracy of the mass accretion rate of the circumbinary disk. 



2.1. Standard gravitational wave phase evolution 



Let us first briefly review the evolution of the GW phase of a coalescing binary, 
neglecting the interaction between the binary and the circumbinary disk. Close enough 
binaries are expected to evolve their orbits by losing the energy and angular momentum 
due to the GW emission. Here, we assume that binary black holes are in a circu l ar or bit, 
because the orbital eccentricity decreases rapidly in proportion to a^^^^^ (jPetersI Il964l ). 
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where a is the semi-major axis of the binary. At the leading order of the Post-Newtonian 
(PN) expansion, the binding energy and orbital angular momentum of a circular binary are, 
respectively, given by 

S = -^^x, J = ^x, (1) 

where c is the speed of light, /i = M1M2/M is the reduced mass with M = Mi + M2 
representing the total mass of binary black holes, x = (GMfiorb/c^)^''^ is a non-dimensional 
quantity of 0((f/c)^), and ilorh = -\/ GM/ is an orbital frequency of the binary. The 
leading terms in the GW energy and angular momentum fluxes emitted from a binary are 
given by the well-known quadrupole formulas as 

• 32c 25 32c 2 5 / n\ 

^gw = -^V X , Jg„ = ^^^^ ^ (^) 

with f] = fi/M and a dot represents differentiation with respect to time. For an equal-mass 
binary, 77 is set to be 0.25. The orbit of the binary gradually decays as a result of the loss 
of its energy and angular momentum due to the GW emission. Equating the time variation 
of the orbital energy or angular momentum with the loss due to the GW emission, one can 
derive, to the lowest PN order, x = GAc^rix^ /5GM, which is translated into the evolution 
equation for the GW frequency / = floj-^/ii as 

/^^ff^V". (3) 



5 

where we introduced the chirp mass. 

Integrating this equation with the coalescence time set to to, the time evolution of the 
GW frequency is computed as 

■T] I 

3/8 



•^^^^ ttGM [5 GM ^ 



9.3xlO-[Hz]r^-2f--'/M^) > (4) 
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where we use 770.25 = ?7/0.25 and rriiQ = M/10Mq for brevity. One can translate this 
frequency evolution into the evolution of a, using the relation (tt/)^ = GM/a^, that is, 

<i) = —ri-. — i^v) (5) 



GM 



tn-t 



1/4 



~ 2.5xlOncm]r,o^/tm,t^^J . (6) 

As the orbit shrinks with time, the GW frequency gradually increases and finally reaches 
the frequency at the inner-most stable circular orbit (ISCO) 

/isco = ^575^ ~ 440[Hz] mro\ (7) 

which is here approximated by the expression in the EMRI limit irrespective of the mass 
ratio. 

The number of the orbital cycles before the coalescence is then calculated as 



to 



Ncyc = J ^dt 

^ 2.3xlO%ol5/W^(^)'''. (8) 

Our main interest is in binary black holes within the mass range 10 — IO^Mq. When the 
mass of each black hole is IOMq and the observation starts one year before the coalescence, 
the initial GW frequency and the initial semi-major axis are ~ O.lHz and ~ 0.047?©, 
respectively. In the typical case one can observe N^yc ~ 10^ cycles of GWs. Such a large 
number of orbital cycles allow us to precisely measure the GW phase correction possibly to 
the level sensitive to the tiny effect of the circumbinary disk. 

The observed Avavefo rm of GWs from the binary at the leading PN order is given by 



(ICulter fc Flanagan 



19941 ) 
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where Q is a factor of at most order unity that depends on the direction and orientation of 
the binary and Dl is the distance from the earth to the GW source. The GW phase (f){t) is 
defined by 

=0o + 27r / fdt, (10) 

J to 

where 0o is the GW phase at the coalescence. 



2.2. Correction to the GW phase due to binary-disk interaction 

In this subsection, we consider the correction to the phase evolution of GWs emitted 
from an inspiralling binary caused by the interaction with the circumbinary disk. Such 
binary-disk interaction induces additional extraction of the orbital angular momentum 
from the binary. This effect will be tiny compared to the phase evolution driven by the 
GW emission when the binary separation is close enough to be observable by future GW 
detectors. Even in that case, however, the correction might be observable because of a 
large number of cycles of GWs in the frequency range of observations. For simplicity, we 
assume that the orbital evolution is dominated by the radiation reaction due to the GW 
emission. In this case it would be natural to assume that the orbit is well approximated 
by a quasi-circular one. Although the circumbinary disk is not necessarily aligned to the 
binary orbital plane, here we also assume a simple aligned configuration. 



The exchange between the orbital angular momentum of the binary, J', and the angular 
momentum of the circumbinary disk is in two ways. One is the angular momentum transfer 



through the tidal/resonant interaction f lArtymowicz fc Lubow 



1994 



Armitage fc Natarajan 



2OO2I ). The other is thro ugh the overflow of the gas from the inner edge of circumbinary disk 



onto the central binary ( lArtymowicz fc Lubow 



1996 



Havasaki et al. 



gas will form accretion disks around respective black holes (iHayasaki et al. 



2007). The transferred 



20081). During 
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this accretion process, some fraction of the angular momentum of the transferred gas will 
temporally go to the accretion disks associated with respective black holes. However, since 
the specific angular momentum of gas at the inner edges of the accretion disks is so tiny, 
only little fraction of the angular momentum of the accretion disks is transferred to the 
black hole spins. Thus, most of this angular momentum is transferred outward via the 
viscous process and i s finally added to the orbital angular momentum of the binary J via 



the tidal interaction flHayasaki 



I2OO9I ). Therefore we can neglect the angular momentum that 
ends up with the spins of the black holes. Thus, the net angular momentum extracted from 
the binary orbital one due to interaction with the circumbinary disk, jTdisk? should be almost 
identical to the net angular momentum flow inside the circumbinary disk, Jdisk, which 
simply follows from the conservation law. (We take the flow of the angular momentum to 
be positive when it flows outward.) 

Now, we roughly estimate the order of magnitude of Jdisk? without relying on a specific 
model of the circumbinary disk. For a quasi-stationary disk, Jdisk should be independent of 
the radius r, which should be equal to the sum of the flows of angular momentum by the 
viscous torque Jvis (> 0) and by the mass accretion Jacc (< 0), 

Jdisk = jvxJyr) + jacc(r) = const. , (11) 

except near the inner edge of the disk located at r = rjn, where the tidal/resonant torque 
is non- negligible. Note that Jacc(^) is negative when the mass is accreting inward. Since 
I Jacc('")| increases with r like oc ^/r for a Keplerian disk, Jvis('") and |Jacc('")| approximately 
balance with each other except near the inner edge of the circumbinary disk. If we were 
discussing an accretion disk around a single central object, the radius of the inner edge of 
the disk would be much smaller than that of the circumbinary disk. In this case, Jvis('") 
should be almost balanced with |Jacc(^)| near r = r^n. In the present situation, however, 
the binary exerts the tidal/resonant torque, which truncates the disk at r = Then, the 
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surface density of the disk is enhanced near the inner edge by a factor of two or so compare d 



Hayasaki et al 



20071 ) 



with the case without the tidal/resonant torque (e.g., see Figure 3 of 
Accordingly, the viscous torque is also locally enhanced by about the same ratio. This 
implies that Jdisk ~ |^acc('"in)|- Thus, the rate of the loss of the binary orbital angular 
momentum is estimated as 



J. 



disk 



J ace ( ^ii 



V2— 



(12) 



where f2in is the Keplerian frequency of the circumbinary disk measured at rin and M is 
set to be positive wh en the mass flux is inward. In the last approximate equality r^n was 



approximated by 2a ( lArtymowicz fc Lubow 



19941 ). Since currently we are not able to make 



the estimate of jTdisk more precise, we introduce a parameter 



3 



Jdisk / M 



J 



(13) 



that represents the efficiency of angular momentum transfer from the binary to the 
circumbinary disk (recall that = rjM). 

The angular momentum transfer from a circular equal-mass binary to the circumbinary 
disk has been recently studied by p erforming hydrodynamic simulations with a-prescription 



( MacFadven fc Milosaylievic 



(jShi and Krolik et al 



2011 



20081 ) and magneto-hydrodynamics (MHD) simulations 



. In both simulations the sound velocity was set by hand 
instead of solving the energy equation. In such a numerical setup, one can freely scale 
the bina ry mass M, the semi-major axis a , and the surface density of the circumbinary 



disk S. iMacFadyen fc Milosavljevid (j2008l ) derived the averaged mass accretion rate 
(M) = 2.5 X 10~^(G'Ma)^/^Eo and the angular momentum flux carried by the circumbinary 
disk (Jdisk) — 1.4 X lO^^GMaSo from the simulated result after 4000 binary orbital 
cycles, wher e Sn is a typical value of S introduced for the normalization. On the 



other hand. 



Shi and Krolik et al. 



mm derived (M) = 1.8 x lQ~'^{GMafl'^T.o and 
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(^disk) — 1-2 X 10~^GMaEo after ^ 77 binary orbital cycles. Substituting these values 
of (M) and (Jdisk) into equation f[T^ . we can estimate the efficiency parameter j to be 
~ 3.6 for the former simulation and ~ 0.69 for the latter simulation, independent of a, M, 
and Sq. The scatter of j between these two simulations will be mainly attribu ted to the 



difference in the strength of shear viscosity. 



MacFadyen fc Milosavljevid (120081) assumed 



that the Shakura-Sunyaev viscosity parar neter ctss to be 0.01, wherea s the effective value of 



ttss derived from the MHD simulation by 



Shi and Krolik et al. 



mm was ~ 0.2. Namely, j 



becomes larger for smaller viscosity. One can explain this tendency qualitatively as follows: 
when the mass accretion rate is fixed, the circumbinary disk becomes more massive for 
smaller viscosity, while a more massive disk gains larger resonant torque from the binary. 

The black hole growth due to the mass accretion from the circumbinary disk to each of 
the black holes via accretion disks can also modify the binary orbital evolution. We note 
from equation f lT3|l that the growth rate of binary mass M/M is the same order as the rate 
of the angular momentum transfer J^disk/<J, if J is the order of unity. Hence, for simplicity, 
we absorb this mass growth effect into the coefficient j in equation ( 1T3|) and neglect the 
time variation of the binary mass in the following discussion. 



From the angular momentum conservation 

~(<-^w ~l~ i^isk) 



dt 



(14) 



we can derive the increasing rate of the GW frequency from equations ([I]), (E]), f lT3|) . and 
(Ml as 



/ 



967r „2 fnGMfY^^ 



-r 



X 



5 J GM _ 



(15) 



32r/2 c3 

The first term in the square brackets of equation f|T5|) is the leading contribution from the 
quadrupole GW radiation given in equation ([3]). The second term represents the effect of 
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the circumbinary disk, and the ellipsis stands for the higher order PN corrections. Recalling 
that X = {tiGM f / c^)^^^ = 0{{v/cY), we note that the correction due to the binary-disk 
interaction has "-4PN" frequency dependence relative to the leading term for a constant j. 

By integrating equation (fT5|) with respect to time, we can express / as a function of t, 
which in turn leads to the expression of t as a function of /: 



5_1 ( 'KGMf \ 



-5/3 



2567r/ V c3 y 
, 5 J GM _4 

1 ^ T^X + 



64 -rf 

The GW phase defined by equation ffTOj) then becomes 



(16) 



16 



TT 



-5/3 



X 



25 J GM _4 

1 7— r — — + 



416 

From the Fourier transformation of equation ([9]), we finally find that the sky-averaged 



17) 



GW waveform in t 



le Fourier domain und er the stationary phase approximation 



(ICulter fc Flanagan 



1994 : 



BertietaL 



20051 ) is given by 



Kf) 



fGM\ 



5/6 



2yT07r2/3 Dl\ J 
where the sky-averaged GW phase: 

71 3 



-7/6gi<E'(/) 



*(/) 



2vr/to 



4 12^ 



TT 



-5/3 



25 J GM 
832^^ 



-X 



(19) 



2.3. Determination accuracy of mass accretion rates: rough estimate 

In this section, we analytically estimate the determination accuracy of the mass 
accretion rate of the circumbinary disk with DECIGO/BBO. The determination accuracy 
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of the mass accretion rate is mainly governed by the low frequency side in the observation 
band, where the effect of the circumbinary disk is the largest. The lower cutoff frequency is 
determined by the duration of the observation or the detector sensitivity, 



/min = max(/(t = tini),/l 



low J 



(20) 



where f(t = tjni) is given by equation (jlj) into which we substitute an appropriate 
observation period to ~ ^ini- We also adopt /low = 10~'^Hz for the cutoff due to the detector 
sensitivity of DECIGO/BBO. Then, /min turns out to be determined by f{t = tini) for 
^0 — ^ini = 5 yr in the mass range 6Mq < M < 5 x IO^Mq. 

The square of the signal-to- noise ratio is defined by 



fmin 



where 



^n(/) 



(21) 



(22) 



and Ndp.v = 8 denot es the effective number of interferometers for DECIGO/BBO 



flCutler fc Holz 



equation (36) of 



20091) and S , 



Yagi et al. 



( f) r epresents a total noise spectral density. As given in 



( I2OIII ). we adopt the non sky-averaged instrumental noise 



spectral density for DECIGO/BBO as 

^;-t(/) = 1.8 X 10"^^ [Hz^^] (^^^ + 2.9 X 10"^^ [Hz-^] 

+9.2 X 10"^2 [Hz"i] 



IHz 

This is obtained by multiplying the s ky-averaged noise spectrum, which is shown in 



(23) 



equation (15) of 



Cutler fc Hod f l2009h . by 3/20 flBerti et al 



2OO5I ). The third term 



of the above equation dominates over the other terms at low frequencies so that 
Sr'if) ~ 9.2 X 10-52 (//lHz)"^ For simplicity, we assume in this subsection that 

/ 



Snif) = 9.2 X 10-52 p^-l] 



IHz 



(24) 



- 14 - 



by neglecting the white dwarf confusion noise from the total noise spectral density. Note 
that the value of /max is irrelevant in the following discussion of this subsection. 

To assess how accurately we can constrain the mass accretion rate, we extract the 
phase correction coming from the interaction between the binary and the circumbinary disk 
in equation (ITUi) . i.e., 

3_2^J_GM_ f^GMfY^^'^ 



*disk(/) - -i^^;^^ [—^ ) ■ (25) 

If there are no degeneracies between the mass accretion rate and the other binary 
parameters, we can roughly claim that the effect of the circumbinary disk is detectable at 
the first order in the mass accretion rate as long as the inequality 

If?(/)W/)I>1 (26) 

is satisfied. This inequality is most likely satisfied at / = /min because \I/disk(/) is proportional 
to a large negative power of /. Therefore, we apply / = /mm to equation (!26l) . As a typical 
source of DECIGO/BBO, we obtain ^(/min) ^ 108 {Dl/3 Gpcy\{to - tim)/5 yr)-i/2 fo^ five 
year observation. Since g{f) is found to be proportional to M^/^/^/^ from equations f|T8|) . 
f l22|) . and f lM|) . it exhibits no dependence on M as long as we adopt /mm = fiUm) oc M~^/®. 
Then, the condition (^E^ gives an analytic estimate for the determination accuracy of the 
mass accretion rate, AM, as 

4 -13/8 5/8 fto — tmi\ '^^^ , , 

4.5 X 10 ""j 'eo.1%25 ^10 (27) 



MEdd ' " '^-^^ V 5yr 

in the mass range 6Mq < M < 3 x IO^Mq, where we normalized the determination accuracy 

by the Eddington accretion rate: 

MEdd = ^ 1.4 X 10'\s-'] e,\m,o , (28) 

and LEdd = 4:nGMmpc/aT is the Eddington luminosity with nip and (Tt denoting the proton 
mass and Thomson scattering cross section, respectively. The mass-to-energy conversion 
efficiency e is set to 0.1 in the following discussion. 
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The relative accuracy in comparison with the Eddington rate decreases with the 
total mass of the binary but increases in proportion to j. This can be understood by the 
facts that the Eddington rate is proportional to M and that ^I/disk(/min) jM~^^^^ in 
equation ( 125|) . which gives AM oc j~^M^^/^. We will compare this analytic estimate with a 
more precise numerical one in the succeeding subsection. 



2.4. Determination accuracy of mass accretion rates: more precise estimate 



We derive the determi nation accuracy of the rn ass accretion rate by applying the 



matched filtering analysis ( ICulter fc Flanagan 



19941 ) ■ We adopt the phase correction 



given in equation ( 1191) as templates, although we do not know the exact form of the 
correction to the GW phase. More detailed study for more reliable prediction of the phase 
correction is left for future work. By neglecting the black hole spins, the determination 
accuracy of the binary parameters 6*" = {In Ai,\nri,tQ,(j)Q, Dl, M} can be estimated by 
Ae- = V(r-^)a„/iVdcv, where 

-/max q}^* I 



ab 



4Re 



is the Fisher matrix. Here, 



fn 



min(/isco, /high) 



df 



(29) 



(30) 



is the cutoff frequency on the higher frequency side, where the cutoff frequency determined 
by the detector sensitivity of DECIGO/BBO is taken as /high = lOOHz and /isco is the 
GW frequency at the ISCO given in equation ([7]). F or the noise spect rum S'„(/), we adopt 



the total noise spectral density (see equation (36) of 



Yagi et al. 



( 20111 )). which includes the 



confusion noise of WD /W D binaries that masks the instrumental noise below / < 0.2 Hz 



( iFarmer fc Phinney 



20031 ). We also assumed that the NS/NS foreground noise can be 



cleaned down to the level below the instrumental noise. 
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Figure [T] shows to what extent one can constrain the mass accretion rate of coalescing 
binary black holes from the GW observation with DECIGO/BBO. The vertical axis shows 
the standard error of the mass accretion rate normalized by the Eddington rate, whereas 
the horizontal axis shows the total mass of binary black holes. The mass accretion rate is 
measurable at a statistically significant level, if it is greater than the determination accuracy 
AM. 

The solid, dashed, and dotted curves represent AM/MEdd for the cases of the 
observation time of Iq — tini =1, 3, and Syr, respectively. The condition AM < Mgdd 
is realized in the mass range 6Mq < M < 3 x 1{]^Mq for one year observation and in 
the mass range 6M0 < M < 8 x lO^M© for both three year and five year observations. 
The determination accuracy is higher for binaries at a closer distance from the earth and 
is also higher with a longer observation time. We also note that the dependence of the 
determination accuracy on the observation time is in good agreement with (to — ^ini)"^^^^ 
shown in equation fl27|) in the mass range 50Mq ^ M < 3 x lO^M©. From Figured! we can 
estimate the mass dependence of the determination accuracy in five year observation as 

^ ^ 1.0 X 10-^j-Vr (31) 

MEdd 

in the mass range 6Mq ^ M < 3 x 10'^ M0. By comparing this estimate with equation (127|) . 
we find that the analytic estimate tends to give about one order of magnitude higher 
determination accuracy and that its dependence on the binary mass is slightly different. 
These are because we have neglected the contribution from the white dwarf confusion noise 
in the analytic estimate. 
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3. Expected number of merger events 

In this section, we estimate the expected number of merger events that have detectable 
signatures of the binary-disk interaction in GWs. 



3.1. Circumbinary disk formation 

The presence of a circumbinary disk around binary black holes with sufficiently large 
accretion rate will not be ubiquitous. To achieve such a large accretion rate as shown in 
equation f lHT]) . it will be required that binary black holes are embedded in a dense gaseous 
environment like a molecular cloud which exists in a star forming region of a galaxy. Its 
typical length scale of the molecular cloud is 10^~^ pc and its molecular hydrogen number 
density is distributed over the range of 10'^~^cm~^. 

If we consider a binary traveling through a dense molecular cloud, the circumbinary 
disk will be formed around the binary through the Bondi-Hoyle-Lyttleton (BHL) accretion. 
The gas within the BHL radius expressed by 

2GM 

tehl = ~ 2.2 X 10-^ [pc] v^Q^m.o (32) 



will accrete onto the binary, where V20 = foo/(20kms~"'^) and foo = y ^'buik + "^s.oo 
with the sound velocity Cs,oo = 0.2kms~^ for typical molecular clouds. Here, we also 
adopted t'buik = 20kms^^ as a fiducial bulk velocity of the binary relative to the ambient 
medium. This is indicated from the observations that the radial velocities of some 
X-ray binaries including GRO J042 2+32 in our galaxy are ~ lOkms"^ (see Table 1 of 



Nelemans fc Van den Heuvell (1l999l )). Since the BHL radius is much less than a typical size 
of a molecular cloud core ~ 1 pc, the circumbinary disk can be formed even in the core of a 
molecular cloud. 
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While the outer edge radius of the circumbinary disk should be less than the BHL 



radiu s, its inner edge radius is approximately given by Tin ~ 2a(t) (lArtymowicz fc Lubow 



1994I ). The consistency condition tbrl > ''^in requires a{t) < cbkl ~ 1-1 x 10~^ [pc] f2^^mio, 
where a{t) can be read from equation ([6]). At five years before the coalescence, this condition 
is easily satisfied. 



The rate of the BHL accretion is given by (iBondi fc Hoyle 



19441 ) 



-3 



~ 2.8 X 10^^ [gs"i] P-i9^^2o'^io , (33) 

where p-ig = poo/(10~^^ gcm~^) with the density of the molecular cloud poo- Its fiducial 
value 10~^^gcm~'^ corresponds to a typical density of dense molecular cloud cores. In 
Figure 1, the dash-dotted line represents this BHL rate normalized by the Eddington rate. 
The determination accuracy is higher than Mbhl in the mass range M > lOOM© for one 
year observation and in the almost whole mass range for both three year and five year 
observations. It can happen that the mass accretion rate is in the range between AM and 
Mbhl; and then the effect of the circumbinary disk is detectable at a statistically significant 
level with DECIGO/BBO. 

From equations fl28p and (155]) . the lowest number density of the molecular cloud for a 
BHL rate to be measured with the determination accuracy AM is given by 

nMc ~ 1-5 X 10^ [cm-3] e^\vl^m^^ { ] , (34) 



Ml 



Edd 



where we substitute AM into Mbhl- 
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3.2. Estimate of merger events 



Following iMii fc Totanil f l2005l ). the volume fraction of molecular clouds with n > umc 



in our galaxy can be modeled as 

(/3-2)(Smc) 



/] 



MC 



5.2 X 10 "e, 



n 



P"'min-'-'g Jn-MC 

10 9/5 -27/5 9/5 



dn 



AM 



-9/5 



0.1 "^20 



10 



(35) 



MEdd 

where n is the molecular hydrogen number density, rimin = lO^cm"'^ and nmax = lO^cm"^ 
are the minimum and maximum number densities of typical molecular clouds, /3 is the 
index of power-law distribution, (Smc) is the mean surface mass density of molecular 



clouds, and is the scale-heig . 
(3 = 2 .8 ( lAgol fc Kamionkowski 



it of the Galactic disk. Here, we adopted t he values of 
2OO2I), (Smc) = 29Mo/pc^ and Hg = 75 pc jsanders et al. 



19841 ). From the condition that nmax should not exceed nuc, we obtain the upper limit of 
AM/MEdd as 



AM 



Ml 



< 6.7 X 10 ^eo.imioW2o^ 



(36) 



Edd 



The merger rate of binary stellar-mass black holes due to the GW emission in our 
galaxy is given by /bh = 2.56 x lO^^yr"^ by "Model A" of Table 4 of 



Belczynski et al. 



(120021 ). The expected number of merger events for the binary black holes within 3 Gpc for 



five year observation is then estimated as 

in 



N 



MC 



D, 



l + z 



to — ^in 
l + Z 



O.U X iU J ' i V20 ?^gal,0.0l"^10 , 



(37) 



where ngai,o.oi = ^^-gai/ (0.01 Mpc~^) with the number density of galaxies ngai and the redshift 
z is set to 0.5 corresponding to Dl = 3 Gpc. We also estimate /mc in the above equation by 
substituting equation f l3T]) into equation (1351) . Note that AM/MEdd satisfies the condition 
given by equation ( l36l) . We have not taken into account the redshift evolution of the 
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merging rate of the binaries per galaxy in the above analysis. When this effect is added, the 
number of events is expected to be a few times larger. We conclude that the feasibility of 
detection of circumbinary disks around binary stellar-mass black holes is marginal in five 
year observation with DECIGO / BBO . We emphasize that our analysis so far does not rely 
on specific models of the circumbinary disk structure. 



4. Astrophysical implications of binary black holes with circumbinary disks 

In this section, we further discuss some astrophysical implications of binary black holes 
with a circumbinary disk. 



4.1. Evolution of binary black holes in the molecular cloud 



The binary transfers its orbital energy and angular momentum to the circumbinary 
disk through the tidal/resonant interaction, which leads to the rapid orbital decay of the 
binary. If the binary evolution towards the coalescence is accelerated significantly during 



its stay in the molecular cloud, the merger rate of binary b 



Belczynski et al. 



ack holes would be increased in 



( I2OO2I ). In this subsection, we 



comparison with the previous estimates of 
examine this possibility. 

Let us consider a circular binary whose coalescence timescale due to the GW emission 



given by ( jPeters 



19641 ) 



a 

64 Vrv 



rgl 
c ?7 



(38) 



exceeds the Hubble time. The binary-disk interaction will be one of the key mechanisms 
to force such binary to merge within the Hubble time. The orbital decay timescale due to 



- 21 - 



the tidal/resonant interaction corresponds to the inverse transfer rate of the orbital angular 
momentum, which is given by equation f[T^ as 

^tide/res = ~ 1-1 X 10^ H J~^eo.l?7o.25"^^^ (39) 



This is consistent with equation (17) of iHayasaki et al.l ( l2010l ) when j = 1.0 is adopted. 



Owing to this tidal/resonant interaction, the orbital evolution is accelerated even when 
tgw > ^tide/res holds, i.e., in the range Otide/res <a< obkl with 

Otidc/res _ / 64 Cttide/rcs \ 

rg V 5 rg 

~ 6.9 X 10^r^/^4/,\;.?5<;/'rh-l/^ (40) 

if the binary stays in the molecular cloud for a sufficiently long term. 

Since the typical size of a giant molecular cloud is tmc ~ 100 pc, the crossing time of 
binary black holes in the molecular cloud is tcross ~ 5 x 10^ [yr] rMc,ioo('ybuik/20kms^^)^\ 
where tmcioo = ''^mc/(100 pc). During this crossing time, the binary can keep interacting 
with the molecular cloud gas. The condition tcross > ^tide/res leads to M > with the 
critical accretion rate 

Mc ~ 2.3 j"^t;2oeo.i'7o.25'"Mc,ioo ^Edd- (41) 

If M > Mc, binary black holes can merge within the Hubble time, even if the initial 
separation is too large to coalesce within the Hubble time by means of the radiation reaction 
due to the GW emission alone. 



4.2. Inner-most disk mass 

In this subsection, we evaluate the mass of the circumbinary disk with the aid of 
the assumption that the structure of the circumbinary disk is described by the standard 
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disk model (jShakura fc Sunyaevlll973l ). The circumbinary disk is assumed to be a steady, 
axisymmetric, geometrically thin, gas-pressure and electron-scattering o pacity dominated 



Keplerian disk. Here, we use the following solution of the standard disk (IKato et al 



20081 ): 



H(r] 



Cs(r) 



4/5 -3/5 1/5 . 3/5 ' 



SS ^0.1 ^10 



4.1 X 10^ [gem' 

1.6 X 10^ [cm] 

-1/10 -1/5 9/10 ■1/5 1 ^ 

1.7 X 10^ [cms^^] 
xa 



-3/5 



21/20 



l/10^-l/5_-l/10 . 1/5 ( J_ 



eo.i "^10 ^ 



-9/20 



(42) 



(43) 



(44) 



where E, if, Cg, rg = GM/(?, and m = M/MEdd are the surface density, the scale- height 
of the disk, the sound velocity measured at the mid-plane temperature of the disk, the 
gravitational radius, and the normalized mass accretion rate, respectively (recall that ass is 
the Shakura-Sunyaev viscosity parameter). 

As it is not the whole mass of the circumbinary disk but the mass at around the inner 
edge of the circumbinary disk that affects the GW phase evolution, we estimate 



M, 



d,in 



7rri„S(rin) 



(45) 



and refers to it as an inner-most disk mass. It can be evaluated by substituting = 2a 
into the expression of S given in equation fH2|) . 



Figure |2] shows the inner-most disk mass of the circumbinary disk that is detectable by 
GWs. They are evaluated at the beginning of the observation period of to — ^ini =1, 3, and 
5 yr, respectively. We also adopt ass = 0.1. It shows that the existence of the circumbinary 
disk can be confirmed, even if the disk mass is less than the binary mass by many orders of 
magnitude. This is because we can observe many (A'cyc ~ 10^) cycles of GWs. 
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We can understand the binary mass dependence of Md^m/M in Figure |2] in the 
following way. The semi-major axis for a fixed value of to ~ ^ini is proportional to M^/^ 
from equation ^ and the surface density at the inner edge of the circumbinary disk is 
proportional to M^/^M^/^a~^/^ from equation ( 142|) . When M is given in equation (127|) . we 
find 

^ ^ 2.2 X 10"i2mtr° (46) 

in the mass range QMq < M < 3 x lO^M© for five year observation. This is roughly in 
agreement with the dotted curve in Figured 



4.3. Decoupling from the circumbinary disk 

When tgfj of the binary is shorter than the viscous timescale evaluated at the inner 
edge of the circumbinary disk, the binary is decoupled from the circumbinary disk. After 
the decoupling, no binary-disk interaction occurs so that it is impossible to detect the effect 
of the circumbinary disk by the GW observation. 

The viscous timescale measured at the inner edge of the circumbinary disk is given by 

2 / \ 7/5 

tvis(a) = ^ ^ tvis,o f - , (47) 

where we used the relation for the disk viscosity u = assCsH and tvis,o = ^vis('^g) ~ 
(l/ass)(crg/Cg(rg)) with equations (l43ll and (l44l) . The decoupling radius at which 
tgw(a) = tvis(a) is given by 

ad f 64 ctvis.o ^ ^^^^ 
-V 



rg V 5 rg 

^ 1 1 V in2n,~^/13 5/13 2/13 1/13 .-2/13 //|o^ 
~ i.i X iU ttgg %25eo.l ^10 ^ ■ l48J 

Even when tg^ is shorter than ttidc/rcs, the binary continues to gravitationally interact 
with the circumbinary disk until the binary orbit decays down to a^. In the mass range 
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of QMq < M < 2 X IO^Mq, the condition a > is satisfied even at one year before the 
coalescence. Namely, the binary continues to gravitationally interact with the circumbinary 
disk, as we described, at least at the early stage of the GW phase evolution during the 
observation period. Thus, although we take no account of the decoupling effect in Figure [H 
this neglect has been justified in the mass range M < lO^Af0. For the larger mass 
range, the determination accuracy is deteriorated by taking into account the decoupling 
effect. However, for M > 1.4 x IO^Mq, DECIGO/BBO is completely insensitive to the 
circumbinary disk when the mass accretion rate is smaller than the Eddington rate. 

The BHL accretion rate can exceed the Eddington rate depending on parameters: 
Poo, Vool and M. X-ray observations have detected the bright X-ray sources , such as 



GRS 1915-1-105 in our galaxy, with a luminosity over the Eddington one (cf. 



Done et al. 



20071 ). Such large luminosities can be explained by a supercritical accretion flow onto 
stellar-mass black holes. If the circumbinary disk is in a super-critical state, equation (|18 
does not apply. Since the accretion timescale is much shorter than the viscous timescale 
in the super-critical state, the binary will not be decoupled from the circumbinary disk 



until the coalescence, a. 



( iTanaka fc Menou 



2010|). 



though it is still unclear how such a circumbinary disk evolves 



5. Summary and discussion 

We have derived the minimum mass accretion rate of a circumbinary disk around binary 
stellar /intermediate-mass black holes necessary for the detection by future space-borne GW 
detectors such as DECIGO/BBO. Our main conclusions are summarized as follows: 

1. The circumbinary disk with M ~ 1.4 x 10^^[gs-i] j-HM/10Mq)33/23 with 10% 
mass-to-energy conversion efficiency, where j represents the efficiency of the angular 
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momentum transfer defined by equation (fT3|) . within 3 Gpc from the earth is detectable 
at a statistically significant level with DECIGO/BBO by observing for five years 
before the coalescence. The determination accuracy of the mass accretion rate is 
higher at a closer distance from the earth, and is also higher for a longer observation 
time. 

2. If binary stellar-mass black holes are residing in sufficiently dense molecular clouds, 
the influence of the circumbinary disk will be detectable. The number of merger 
events of such binaries within 3 Gpc for five year observation, is approximately 
estimated to be (9(0.1) by assuming typical values for the event rate of the black hole 
coalescence and that molecular clouds have the number density from 10~^ cm"'^ to 
riMc given in equation (134|) . This suggests that there might be a possibility to detect 
the circumbinary disk around binary stellar-mass black holes with DECIGO/BBO. 



In the estimate presented in this paper, we still have a lot of uncertainties, e.g., the 
formation rate of binary black holes, the molecular hydrogen number density, the relative 
velocity to the molecular cloud, the efficiency of angular momentum transfer (j: see 
equation f|T3|) ). and so on. Therefore, we should understand that there remains a possibility 
that the event rate is one or two orders of magnitude larger than our estimate. 

For the idealized situation that we discussed in this paper, the correction to the 



GW phase evolution shows — 4PN frequency de pendence. 



as the one for a braneworld model discussed in 



acceleration ( jSeto et al. 



Yagi et al 



his d ependence is the same 



20111). Similarly, the cosmic 



200ll ) and the alt ernative theories o f gravity in which the 



gravitational constant G is not a costant flYunes et al 



20101 ) also give the corrections 



with — 4PN frequency dependence. Furthermore, the acceleration acting on a binary in 
time-indepe ndent external gravi tational field gives the correction with the same frequency 



dependence (lYunes et al. 



2011al ). Therefore, once the deviation from the standard template 
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is discovered, discriminating the effect of the circumbinary disk from the others will be a 
difficult task. However, for example, the binary-disk interaction rapidly decays after the 
decoupling of the circumbinary disk as we discussed in section 14.31 This is a signature 
unique to the interaction with the circumbinary disk. In addition, there can be a method 
to distinguish the effect of the circumbinary disk from that of modified gravity theory by 
detecting the GW signals from multiple sources. This is because the effect of modified 
gravity theory is universal, whereas that of circumbinary disk is not. Therefore, further 
detailed studies may provide some ways of distinctions. 

In this paper we have assumed that the binary is in a circular orbit. The orbital 
eccentricity also gives the frequency dependence of negative PN order (-19/6PN) 



(ICutler &: Harms 



20061 ). Inclusion of the orbital eccentricity as a parameter of the waveform 
reduces the determination accuracy of — 4PN order term by some factors. However, we 
expect that we can determine both the orbital eccentricity and the mass accretion rate 
for the large rn a ss par ameter range with errors of (9(0.1) for typical binary parameters. 



Key fc Cornish! (120111 ) calculated the determination accuracies of full seventeen binary 
parameters including an initial orbital eccentricity for binary massive black holes with 
LISA, using a Bayesian analysis under Markov Chain Monte Carlo simulations. They 
conclude that the LISA could distinguish between circular orbit and eccentric orbit with a 
very small orbital eccentricity e ~ 10"'^. 

Furthermore, we would like to point out that the effect of the circumbinary disk 
can be important, even if it is suppressed to the undetectable level for each event. The 
primary mission of the space-borne GW antennas such as DECIGO/BBO is to measure the 
primordial GW background radiation. In order to achieve this goal, abundant binary GW 
sources become the foreground noises. To remove them, the expected GW signals for the 
events identified by the matched filtering will be subtracted from the data. However, if the 
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theoretical templates contain some unkno wn elements, the subt r action leaves residu als that 



remain to contribute to confusion noises (jCutler fc Harms 



2006 



Harms et al. 



20081 ) ■ 



Finally, we briefly mention the relevance of our discussion to the ground-based GW 
detectors. With these detectors, the distance to which we can detect GWs from binary 
black holes is reduced, and hence the expected event rate is extremely small. Using the 
third-generation- ground-based GW detector: Einstein Telescope (ET), and assuming high 
signal-to- noise ratio corresponding to an accidentally nearby event at -D^ = lOOMpc, the 
effect of the circumbinary disk will be detectable only for the mass accretion rate 10"^ 
times higher than the Eddington rate. Even if stellar-mass binary black holes are at 
rest in a dense molecular cloud core, the expected BHL accretion rate is 10'^"^ times as 
much as the Eddington rate. Hence, we would be able to safely conclude that the effect 
of the circumbinary disk becomes relevant only for space-borne GW detectors such as 
DECIGO/BBO. 
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Fig. 1. — Determination accuracy AM of mass accretion rates of the circumbinary disk 
around coalescing binary black holes within Dl = 3 Gpc. Here, j = 1 is adopted (see 
equation ( fT3l) for the definition of j). The determination accuracy is normalized by the 
Eddington accretion rate with 10% mass-to-energy conversion efficiency. M is the total 
mass of binary black holes in units of the solar mass. The solid, dashed and dotted curves 
correspond to the determination accuracies with the observation time, to — tmi =1, 3, and 
5yr, respectively. The dash-dotted line shows the Bondi-Hoyle-Lyttelton accretion rate 
normalized by the Eddington rate, where the BHL rate is given by equation ( 133|) . The mass 
accretion rate can be measured from GW observations at a statistically significant level if 
AM is less than the fiducial accretion rate M. The horizontal solid thin line represents 



10-12 I I 

1 10 100 1000 10* 10° 

Fig. 2. — The inner-most disk mass that gives detectable effect on the GW phase evolution 
is plotted as a function of the total mass of binary black holes. The vertical axis is normalized 
by the total mass of binary black holes, and the inner-most disk mass is defined by in = 
7rrfnS(rin), where rjn and S are the inner edge radius of circumbinary disk and the surface 
density of the circumbinary disk evaluated at around Tin, respectively (see equation fjlSj) in 
the text). The solid, dashed and dotted curves correspond to M^^^^/M for to — Uni =1, 3, 
and Syr, respectively. 
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